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ABSTRACT 

> : 

All the next-to-leading order contributions to the quasi-particle dispersion 
laws of a quark-gluon plasma which due to infrared singularities are sensitive to 
the magnetic-mass scale are calculated using Braaten-Pisarski resummation. 
These relative-order-g ln(g) corrections are shown here to generally contribute 
to the dynamical screening of gluonic fields with frequencies below the plasma 
frequency as well as to the damping of propagating gluonic and fermionic 
quasi-particles. In the limit of vanishing wave- vector the infrared singularities 
disappear, but in a way that raises the possibility for formally higher orders 
of the Braaten-Pisarski scheme to equally contribute at next-to-leading order 
when the wave-vector is of the order of or less than the magnetic-mass scale. 
This is argued to be a problem in particular for the fermionic damping rate. 
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1. Introduction 



The leading-order results for the dispersion laws of quasi-particles in a quark-gluon 
plasma |l] are well known and are readily obtained from the high-temperature limit of one- 
loop Green's functions or from solving the collisionless Boltzmann equation [H. However, 
anything beyond the leading terms becomes accessible only through a resummation of the 
conventional perturbation series. 

A systematically improved perturbation theory, which turns out to involve single pow- 
ers of the QCD coupling constant g rather than the usual g 2 , has been established some 
years ago by Braaten and Pisarski ||. It is based on a resummation of all the leading-order 
self-energies and vertex functions furnished by the so-called hard thermal loops @, H] • In 
a first application this led to the solution of the long-standing problem H of how to cal- 
culate the damping constant of the lowest QCD plasmon excitation with vanishing wave 
vector in a gauge- invariant way ||. In later work it was found that with non-zero wave 
vector the damping of gluonic as well as fermionic excitations is infrared divergent with a 
logarithmic singularity which can be cut off by a finite screening mass for static magnetic 
fields. This gives rise to a contribution proportional to ln(m e i./m magn .) ~ ln(l/g) whose 
coefficient is calculable perturbatively |7|, ||, [|, |I0|, Similar singularities have also been 
found recently in the next-to-leading order calculation of the nonabelian Debye screening 
mass 0, |13|] from the pole of the static gluon propagator (at imaginary wave vectors) as 
well as in the perturbative evaluation of the correlator of Polyakov loops |13|, |T^|| . 

Such a sensitivity to the magnetic mass scale comes as a surprise since by superficial 
infrared power counting |15|] one would not expect it already at (resummed) one-loop 



order. It is in fact due to the necessity to evaluate the loop diagrams at the position of 
the pole of the leading-order propagator, which leads to "mass-shell" singularities in the 
presence of the massless magnetostatic mode. Indeed, these singularities appear also in 
the case of QED, where no magnetic screening mass can be generated. Presumably the 
finite width 7 of the full propagators is also able to provide the required infrared cutoff 
H, which would again lead to ln(m e i./7) ~ ln(l/g), or ln(l/e) in the case of QED. 

In this work we shall present a unified treatment of all these next-to-leading order 
corrections to gluonic and fermionic dispersion laws which are infrared singular due to 
mass-shell singularities (Section 2) and therefore sensitive to the 'magnetic-mass scale', 
by which we mean simply the scale of new physics that acts as an effective infrared cutoff, 
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be it an actual magnetic mass or the shielding of the singularities through a finite width 
of the full propagators or another mechanism. For moving quasi-particle excitations we 
reproduce the results of Pisarski |TTJ who found that the infrared-singular contributions 
are proportional to the group velocity of the respective modes. For the case of dynamical 
screening of perturbations with frequencies below the plasma frequency we find direct 
proportionality to the modulus of the (now imaginary) wave vector (Sections 3 and 4). 
We also take a somewhat closer look (Section 5) at the limit of vanishing wave-vector 
where the infrared singularities seem to be absent and we shall argue that for wave- 
vectors of the order of or smaller than the magnetic-mass scale these infrared singularities 
could still leave non-negligible imprints which would render the resummation scheme of 
Braaten and Pisarski incomplete in this limit. 



2. Quasi-particle mass-shell singularities 



The leading-order finite-temperature corrections to the gluon and fermion self-energies 
give rise to effective (albeit momentum-dependent) thermal masses of the order gT, 
where g is the QCD coupling constant and T the temperature. Because of a nontrivial 
tensor structure, they give rise to different dispersion laws for spatially longitudinal and 
transverse gluons according to 

Q 2 = n e (Q ,q) 
Q 2 = n t (Q ,q) 



\ (3- 2 - n t ) 



where IJ^ is the gluon self-energy and m 2 = (N + Nf /2)(gT/3) 2 for gauge group SU(iV) 
with Nf flavors. For fermions (whose zero-temperature rest-mass we assume to be -C gT), 
there are also two different modes at finite temperature from the solution of 

det(#-E) = , (3) 

where £ is the fermion self-energy. Eq. (F|) has two solutions for positive q, which are 
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Figure 1 : The leading-order gluonic and fermionic dispersion curves and their infrared 
singularities at subleading order. For the transverse (t) and longitudinal (£) 
dispersion curves of the gluonic modes UJ and q are given in units of m ; for 
the quark (+) and plasmino (— ) modes the unit is M 2 . 

given by 

Qo = ±q + n ± (Q ,q) (4) 

with n^^\(^9AjQi±l) ±2 } (5) 

2g [\ g J \Qo-qJ \ K J 

and M 2 = C F {gT) 2 /8, C F = (N 2 - 1)/(2N). 

In Fig. 1, the dispersion laws Qo = ou(q) of the above modes are shown as curves in 
the uj 2 , g 2 -plane. For real q, q 2 > 0, they correspond to propagating quasi-particles; for 
q 2 < 0, they describe the screening of fields oscillating below the plasma frequency, where 
\q\ is the inverse screening length. 

The physical significance of the various modes have been discussed in full length by 
Weldon JTJ, to which we refer for further detail. Let us just mention here that the spa- 
tially longitudinal gluonic excitation (which is often referred to as the plasmon) and the 
fermionic one corresponding to the lower sign in Eq. @) (sometimes called the plasmino) 
have no counterpart at zero temperature. They can only be understood as collective 
phenomena. 

If one attempts to go beyond the leading-order results just presented, one has to face 
the problem that the conventional, bare perturbation theory becomes insufficient when- 
ever the loop momenta probe the new ( "soft" ) scale ~ m, M. This is true in particular 
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for the next-to-leading order terms in Eqs. (|l]j2|]4]) . An improved perturbation theory, 
due to Braaten and Pisarski ||, requires a resummation of both the leading-order con- 
tributions to the various self-energies and also to the vertices, collectively termed 'hard 
thermal loops'. However, in a perturbative treatment of the solutions to Eqs. (|l],0,fD one 
has to evaluate the correction terms to the self-energies at the location of the poles of the 
leading-order propagators, which in the presence of massless modes generally will give rise 
to "mass-shell" singularities. Indeed, while all the propagating physical modes have non- 
vanishing effective masses, the thermal mass of the transverse gluon propagator vanishes 
in the static limit, corresponding to the absence of screening of static chromomagnetic 
fields up to distances ~ l/(gT). 

Consider the following integral which appears in the resummed one-loop contribution 
to the self-energy of mode i, where in the loop integral the propagator for mode i is 
coupled to the one of transverse gluons, 



<Si(Qo,q) 



E 

P =2mnT 

E 

P =2ninT 



d 3 p 
(2^)3 

d 3 p 



A t (P) MQ-P) 



1 



(2VT)3 P 2 _ n t ( P) ( Q _ p)2 _ n .(Q _ P) 

For gluons, i either means t or £; for fermions the second propagator in ( 

[K -k-n + (K)] [K + k-n_(K)] = Af{K) 

with K = Q- P. 



• (6) 

Q 2 =n t {Q) 

|) has the form 

(7) 



The n = contribution in the above sum comes with a vanishing n t (Po = 0,p) = 0, 
and for small p the second propagator is ~ l/(pq) on the mass-shell of mode i, so 
this contribution involves a logarithmic singularity on mass-shell. Concentrating on the 
infrared-singular part, we therefore have 



d 3 p 1 



1 



(27r) 3 p 2 p 2 — 2pq + 
T [l + d^n^q)}- 1 J 



d 3 p 



P) 



Pii{Q) - iea 
1 



Po=0 



(27r) 3 p 2 p 2 — 2pq — iea 1 



where here and in the following ~ indicates that we have dropped regular terms. In the 
last line of Eq. @ we have used that the difference of the two i7j when taken at Pq = 
is a function of Q and p 2 — 2pq. Only the first term of its Taylor series with respect 
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to the latter variable is relevant for the infrared singular contribution; the others lead to 
regular integrals of the form fd 3 pp~ 2 (p 2 — 2pq) n , n > . In Eq. (H) a is a sign fixed in 
accordance to retarded boundary conditions by Q — > u> + ie: 

a = sign (u [1 -d^II^u.q)]) (9) 

which leads to 

a' = sign (u [ 1 - d^II^u, q) ] / [ 1 + d^II^u, q) ) ) = sign (v*) (10) 

with fj the group velocity dui/dq of mode i. Note that —uP'd^ = q 2 d q 2 since IIi(uj,q) 
depends only on c<j/g. 

In order to give meaning to the singular expression (f|), we introduce an infrared cut-off 
A <C gT for p : 

, ; A) = * [ 1 + W , ,) ]- f°° * i- In ( " + *> - fe , (11) 
47r^ J\ p 2q \p — 2q — tea' J 

which has a regular real part but a singular imaginary part, where A cannot be removed, 
S l (u,q;X)=ta'^[l + d q2 n t ]- 1 ln(^j+O(\ ) . (12) 

For frequencies below the plasma frequency m, where the gluonic dispersion laws 
describe screening corresponding to poles at purely imaginary values of q = ±i\q\, we 
have instead 



Si(uj , q = ±i\q\; X) = [ 1 + cbiJ; ] 1 [ — arctan f — ) 

47r^ 7a p |(?| \ P J 

= 84l |1+a ' 277 ' rl ' n (T) +0 ( A °) (13) 

which now has a singular rea/ part. 



The results ( P^JT3| ) are only valid when |g| ^> A. In the limit |g| — > 0, 5 is in fact 



linearly divergent rather than only logarithmically, which makes S proportional to 1/A: 
Si(u, q; A) -> ^ [ 1 + a 9 2i7i ] _1 + O f ^|) when \q\ < A . (14) 



4tt 2 A l ' 9 tJ ' VA 

In the applications below it will turn out that all the infrared singular contributions 
to the next-to-leading order self-energies involve rational functions of the momenta times 
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the paradigmatic expression S such that the limit q — > is infrared finite and, apparently, 
the linear divergence in fll4] ) is defused. However, we shall find reason to reassess the case 
q — > more cautiously later, so let us for the time being restrict our attention to the case 
of a non-zero q. 

With non-zero q ~ gT, the infrared singularities are only logarithmic and it suffices to 
assume that an actual infrared cutoff will be generated at some scale \x ~ g x T, x > 1 to 
make sense of the calculations. The technical cut-off A can then be identified with fi as far 
as the logarithmic term is concerned, which turns into (x— 1) \n(l/g). We expect /i ~ g 2 T 
and shall refer to it loosely as the magnetic-mass scale, but actually we mean simply the 
mass scale where new (infrared) physics eventually removes the mass-shell singularities of 
Eq. (|8|). This may be through the nonperturbative generation of a screening mass for static 
chromomagnetic fields or something else, e.g. the finite width of the true quasi-particle 
excitations |J, which, as we shall see, is again of the order of g 2 T (up to logarithms ln(g)). 

At any rate, this assumption fixes the term involving ln(l/g), while the finer details of 
the physics pertaining to the supersoft scales, which are outside the scope of this paper, 
will be relevant only for the subleading contributions. 



3. Gluon self-energy 



In resummed perturbation theory, the one-loop correction to the gluon self-energy is 
usually a rather complicated object, since also the vertices have to be dressed. For non- 
zero external wave- vector the resulting algebra is considerable. Some of its properties have 
recently been studied by two of the present authors [[16]. Fortunately, when concentrating 
on the infrared-singular parts, this algebra can be greatly simplified and we shall in the 
following describe only the shortcut to obtaining them. In an appendix we display the 
next-to-leading order self-energy in full detail, from which the following results can be 
straightforwardly reproduced, albeit in a less transparent manner. 

First notice that only the diagrams containing two propagators are able to produce 
mass-shell singularities, of which one has vanishing 4- momentum. The most singular terms 
are those that have no further loop momenta in the numerators, so all the momentum 
algebra factorises into a prefactor to singular integrals of the form In particular, 
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only the part of the complicated dressed vertices contributes where one leg has vanishing 
momentum. Thanks to the gauge invariance of the hard-thermal-loop effective action || , 
these vertices are determined by differential Ward identities, 

*IW-Q, Q, 0) = d Q , *n^{Q) (f 5) 

where the vertex *r is the sum of the bare and the hard-thermal-loop vertex and *Z7 the 
sum of the bare kinetic term and the hard-thermal-loop self-energy : 

* n i»> = qvqv -Q 2 g ^ + = (n t -Q 2 )A^+(n t -Q 2 )B^ . (16) 

Here A is the projector on spatially transverse vectors (s. fllTf ) below) and B^ u = g^ v — 
Q»Q V /Q 2 - A» v . 

The full next-to-leading order gluon self-energy is in general gauge-parameter depen- 
dent and non-transverse, but one can show on an algebraic level that the gauge parameter 
drops out from 5II tt e on the respective mass-shells and also that SII^ is transverse (in 



the four-dimensional sense) on the longitudinal mass-shell [|T(J. This on-shell gauge inde- 
pendence allows us to further simplify the algebra by using Feynman gauge, where indeed 
only the integral (JD) arises, and from on-shell transversality we have 



5n t (Q) = ^A^6n itv = ~\6 mn -^)6n mn (17) 



sn e (Q) = -^sn 00 . (is) 

The whole algebra is now a matter of a few lines only: 



8II t (Q) ~ (s mn - ^ ) *T mai (-Q, Q, 0) T ara (-Q, Q, 0) S t (Q) 



-g 2 N (d.At 1 ) 2 S t (Q) 

-g 2 N4q 2 [ 1 + d?II t (u, q) f S t {Q) (19) 



and 



m,(Q) - -9 2 N[%\ T^(-().0.())cS ; f()) 



q 2 



-g 2 N (d qi Aff S t {Q) 

-g 2 N4q 2 [l + d q2 n £ (uj,q)} 2 S i (Q) . (20) 
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This exhibits the prominent role played by the factors Si, i = t,£, studied in the previous 
section. As they have no singular real parts at positive q 2 and no singular imaginary part at 
q 2 < 0, (0) and (|2D|) tell us that there is no danger from the infrared in the corresponding 
parts of the polarization functions at next-to-leading order. For the singular parts, on the 
other hand, S contains in one simple expression the essentials on both next-to- leading- 
order damping and screening for q 2 > and q 2 < 0, respectively (see Figure 1). 

Let us first consider the consequences for propagating quasi-particle excitations which 
have uj > m and q ~ m 3> \i. Inserting the on-shell values of the derivatives of the various 
iZ; we obtain 

6n t ~ -*9 2 NT^- [ 1 + d q ,n t (u, q) ] In Q 
Sn,(Q) ~ -ig 2 NT^ [ 1 + d#n,(ui, q) ] In 

= -'^ T l (™ 2 - « 2 ) \* g) ■ m 

where in evaluating S we have assumed m/fM ~ 1/g and used that a' = +1 for both 
modes. 

In order to derive corrections to the dispersion laws, one has to take into account 
that the leading-order U t i also varies when uJt,i(q) changes. Expanding the condition 
u> 2 e (q) = q 2 + II tj i + 5 lit/ around the leading-order result u>o(q) 2 , we have 

bu\ t = (uj 2 - U)l) u = 5II t/ (uj , q) + 5uj 2 t/ d u 2lI t/ (uj Q , q) +0 (V) 

= i 7'n'f i + • < 23 > 

I - d^n^Uo, q) V J 

Interpreting the imaginary correction 5uj 2 = — 2iu;o7 as a damping term 7, we obtain 
(dropping the index on uj in the final results, since the difference is of higher order now) 

9 2 NT q [l + d q2 n t/ (u,q)} f l\ 
7M 4vr u[l-d^n t)t (u,q)] * \g J 
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with 

io 3(m 2 - Q 2 ) J lo 3(m 2 u 2 - Q 4 ) 

Vl = ~ q 3m 2 -Q 2 Vt = ~ q 3m 2 u 2 -Q 4 ' (25) 

both of which are zero at q = and approach 1 with increasing q. This exactly reproduces 
and confirms the results of Ref. JTTJ . 



The above calculations allows us equally to derive the correction to the part of the 
dispersion laws which describe screening of external perturbations with frequencies below 
the plasma frequency, ui < m (see Figure 1). There the gluon propagator exhibits poles 
for imaginary values of q, which gives the effective screening length for a given frequency 
to. In the static limit this reduces to the electric (Debye) mass and the magnetic mass 
(which is zero at leading order) for mode I and t, respectively. In the general case, both 
screening lengths are finite (which is sometimes referred to as "dynamical screening"). 
Here it is more natural to keep to fixed and to determine Sq(u) according to 

1 + d^II ti t{u, q) v ' 
with q{to) as given by the leading-order results. 

Together with the results for St/ for imaginary q, this gives for \q\ ~ m ^> \i 

%l,^.n(l) , 

so the correction is now a simple constant. In the static limit \q\i = \f^m = m e \. and 
(p7| ) agrees with the next-to-leading order result for the electric (Debye) screening mass 



of Refs. 12, 13, 14 



In the previous section we have seen that the logarithmic singularity of 
g 2 S ~ (g 2 T/q) ln(q/n) ~ g\n(l/g) for |g| ~ m 
turns into a linear one for |g| ^0: 

g 2 S ~ gm/jj, for \q\ < n , 



which for // ~ g 2 T is 0(1). The final results for the gluon self-energy ([21], 22) come with 
a prefactor that vanishes like q 2 , so that both, the damping coefficient and the correction 
to dynamical screening, as given above, do not receive contributions from the singular 
integral S. This could change, however, when higher-order corrections are taken into 
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account, and we shall resume this point after extending the above calculations to the 
fermionic spectrum. 



4. Quark self-energy 



The calculation of the infrared- singular parts of the one-loop resummed quark self- 
energy can be simplified in full analogy to the gluonic case of the previous section. Again 
the result factorises into an expression involving only the external momentum and the 
singular integral Sf introduced in sect. 2. The dressed quark-quark-gluon vertex is needed 
with vanishing gluon momentum only and can again be derived from the self-energy 
through a differential Ward identity 

*T li (-Q,Q,0)=d QI *($-E(Q)) , (28) 

so that 

= -g 2 C F [d qm ($-Z(Q))} (#-E(Q)) [<($-£(£))] S f . (29) 

The corrections to the two branches of the fermionic dispersion laws (see Figure 1) are 
determined by 

5U ± = itr (( 7 ° =F i/q)SZ) . (30) 

Introducing 

D± = cu T q - n± (31) 

we obtain 

5n ± = - \D T {d q D ± f + l -D ± (D + - D_) 2 } S f , (32) 

where only the first term in the curly brackets survives on the respective mass-shell 
D± = 0. 

In the fermionic case, the propagator has poles only for real values of q, corresponding 
to the fact that fermionic quantum numbers cannot get screened. For q ~ M ^> ji ~ g 2 T 
we have 

9 2 C F T 
4vr '-"-^'-Vo 



§n ± = -i \ \d q D±\ ln( - 



. g 2 c F T 

An 



2ou(Q 2 - M 2 
IT- ' 



qQ 2 



ln(l) , (33) 
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where in evaluating Sf we have used 



1 + d a 2ll 



f 



1 

2q 



--d q (D+D.) 



(34) 



and sign(v-i-) = — sign(d q D±) in evaluating Sf. Determining the correction term to the 
dispersion law uj±(q) = uo±(q) + 8u±(q) requires again to take into account the induced 
variation of the leading-order self-energy through 5u± according to 



sn ± (u ,q) 



6II±(uJo,q) 



1 - II± (uj , q) d Uo D±(u ,q) 
Since <5i7± is purely imaginary, this yields a fermionic damping coefficient 



(35) 



7± 



g 2 C F T \d q D ± \ 



In 



4tt d w D±(uj,q) \g 

9 2 C F T 
4ir 



uj{Q 2 - M 2 ) T qQ 2 



qM 2 



In - 



(36) 



in accordance with the results of Ref. IffTI . 



The group velocity v± equals ±1/3 in the limit q — ► 0, monotoneously increasing 
towards +1 for large q (with a zero for the plasmino mode at q ~ 0.41M). For q — > 0, 
(P^) is no longer valid but becomes proportional to qj \i in place of ln(l/g), so that the 
infrared- singular contribution disappears. This agrees with the calculations performed in 
Ref. |I7| at strictly q — 0, where a finite result without ln(l/g) was obtained. 



5. Higher-order contributions 



For q 0, we had to invoke higher-order corrections for providing a physical infrared 
cutoff for the mass-shell singularities in the next-to-leading order corrections to the var- 
ious on-shell self-energy components. Since the singularities were only logarithmic, the 
assumption that such an infrared cutoff will indeed be put into effect by the higher-order 
contributions at a certain scale /x suffices to fix the dominant term oc gln(l/g). Clearly, 
the other terms of 0(g) depend on the details of the new infrared physics at the scale 
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jj, and thus are beyond the perturbative scheme of Braaten and Pisarski, which is based 
on a resummation of the leading-order terms pertaining to the scale m ~ gT. Only the 
real-part corrections to the dispersion laws at 0(g) are infrared-safe and calculable this 

way 



18, 16 



In the long- wavelength-limit, q — > 0, this sensitivity to scales <C m disappeared because 
the potentially singular terms in the self-energy corrections were of the form 

5n tt i(u,q) ~ q 2 ft,e(u/q)g 2 S t ,£ (37) 
Sn ± (u,q) ~ qf ± (uj/q)g 2 S f . (38) 

In sect. 2 we have observed that S becomes linearly singular for q — > 0, and assuming 
that higher-order corrections render S finite so that S ~ T/u, one can indeed drop the 
above contributions at q — 0. 

However, this may cease to be justified if there are higher-order corrections to the 
prefactor of S that do not vanish when q — > 0. In the resummed perturbation theory, 
one collects all the contributions at the soft scale gT, consistently disregarding potential 
terms in the effective action that are proportional to fi/m and therefore suppressed by 
extra powers of g. So higher-order corrections to the dressed vertices and propagators 
could in principle change the prefactor in (1371) like 



q — > q + Cigmq + c 2 gmuj + c 3 gm + . . . (39) 



and in (|38| ) like 

q — >• q + digm + . . . , (40) 

with dimensionless functions Ci(u,q) and di(u,q) . 

If fi which cuts off the linear singularity of S for q < \i was much smaller than gin, 
then there would even be the possibility that the leading- order results might get modified 
by the higher-loop orders, but only for q < \l. We shall exclude this rather improbable 
eventuality by assuming that \i ~ gm ~ g 2 T ' . 

Let us first consider the effect of nonzero c's on the result for dynamical screening, 
eq. (|26|), for q < \x. The correction to the inverse screening length \q\ becomes 

ii ii i s 

gmq g m g m u g m 

°\q\t,t ~ Hc x hc 2 hc 3 h... 

fj, fj, qfi qfi 

gmu gm 2 

~ gq + Cigm + c 2 h c 3 \- . . . (41) 
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Obviously, 02(0;, q) and cs(u>, q) should vanish for q — * in order that the kinematical 
situation for q — > does not become singular, which we shall take for granted. 

In the static limit, the transverse branch of the dispersion laws has \q\ = at leading- 
order, i.e. a vanishing magnetic mass. A nonvanishing 5\q\ t for uj —>■ would be in- 
terpreted as the generation of a magnetic mass m magn . The result obtained within the 
Braaten-Pisarski scheme vanishes in this limit, which is in agreement with the null result 



of Ref. [|T2|. But a nonzero C\ would render m magn ~ Cigm, which is consistent with our 
assumption that the cut-off n is of the order of the magnetic mass (although not neces- 
sarily identical with it). The linear mass-shell singularity of S for q < 11 could thus play 
a prominent role for the generation of a magnetic screening mass through higher-order 
corrections. 

Turning now to the propagating gluonic modes, we note that for q — > 0, the singular 
contributions to the damping constant still vanish like 

gq 2 g 2 qm 

lt,e ~ h ci h ... , (42) 

since we have excluded C2 and C3 for q — > 0. Thus there seems to be little danger that 
the relative-order-g results that have been obtained previously for q = || [18|] could be 
modified by higher-order corrections. Also for nonzero q < fi, the contribution from the 
mass-shell singularities remains below O(gm). 

The situation is somewhat different for the fermionic modes, however. There we have 

gqm g 2 m 2 

7± < + di + ••• , (43) 

/i 11 

and already the first term, which arises within the Braaten-Pisarski scheme, is of the order 
gm ~ g 2 T for q ~ 11. For such momenta 7 is obviously not calculable within the one-loop 
resummed approximation, because through the linear divergence of S, loop momenta of 
the order of /i contribute on a par with the ones of order m. For q — ► 0, these contributions 
are suppressed, but if higher-order terms could produce a nonvanishing d\, then also the 
strict q = result for 7 at order g 2 T would become infested by higher-order contributions. 

However, the proportionality of 5II± to q can be traced back to the Ward identity 
(p8|) which is responsible for the simple form of (j32|): therein both D + and D_ vanish 
like q for q — > 0. In the Abelian case, the tree-level-like Ward identities hold also beyond 
the level of hard thermal loops, whereas in the nonabelian case, one can retain simple 
Ward identites by choosing to work in axial gauge, so it seems plausible that d\ oc q, 
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which makes the results obtained for strictly q = in Ref. []17] stable against higher-order 
contributions, despite the lurking mass-shell singularities. 



6. Conclusion 



We have found that the next-to-leading order corrections to dynamical screening 
masses for gluonic fields with frequencies below the plasma frequency as well as the 
damping constants for the propagating quasi-particle modes are strongly sensitive to the 
magnetic-mass scale except when the wave-vector is exactly zero. On the other hand, the 
next-to-leading order corrections to the real part of the dispersion laws of the propagating 
modes have turned out to be infrared- safe. 

The infrared singular contributions for both, screening and damping, are determined 
in essence by one simple expression exhibiting (quasi-particle) mass-shell singularities. 
For large enough modulus of the wave-vector the latter are logarithmic, yielding contri- 
butions of relative order g\n(l/g). The coefficient in front of the logarithm is calculable 
once the scale of the cutoff brought about by higher-order contributions has been deter- 
mined. The coefficients under the logarithm ln(l/g), which are of obvious importance 
when the coupling constant is not infinitesimally small, are clearly beyond these pertur- 
bative considerations. 

For small wave-vector, we have found that the mass-shell singularities become even 
linear, which in principle opens a way for higher-order corrections to contribute on a par 
with the ones obtained within the resummed perturbation theory of Braaten and Pisarski. 
We have argued that it is plausible that the case of exactly vanishing wave-vector is 
stable against such corrections, whereas the particular case of the damping constant of 
fermionic modes with nonzero wave- vector of the order of the magnetic mass-scale remains 
uncalculable. 

Linear mass-shell singularities have previously been identified as the root of a potential 
problem with gauge independence of next-to-leading order corrections to the disper- 



sion laws. Whereas formally one can prove gauge fixing independence ||20|| , in covariant 
gauges the unphysical modes of the gluons behave like zero-mass particles, and lead to 
linear divergences in the residues of the quasi-particle propagators. Unless they are iden- 
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tified as such by the introduction of a (purely technical) infrared cut-off, they can mimic 
contributions to the pole position pT| . 



However, the type of mass-shell singularities that we have discussed in this paper 
appear directly in the corrections to the pole position. They call for a physical cut-off to 
be provided by higher-order corrections. Indeed, in the results for dynamical screening of 
transverse gluonic modes, we have seen that the linear mass-shell singularities themselves 
could play an important role in a dynamical generation of a magnetic screening mass 
~ g 2 T, which would be the most obvious candidate for such a physical infrared cut-off. 

In the nonabelian case, the mass-shell singularities are there even in the purely static 
situation, where they provide the dominant next-to-leading order contribution to the 
Debye screening mass. Here the other possibility for an effective infrared cut-off that has 
been discussed in the literature, namely damping of the internal propagators, can hardly 
be operative, which underlines the need of (chromo-)magnetostatic screening. 

On the other hand, in the Abelian case, a magnetic screening mass cannot be gen- 
erated. This is no problem for the Debye mass, which is infrared safe in QED. But 
mass-shell singularities are there in the next-to-leading order corrections to the electron 
propagator, and there has been some controversy [22, |24[ on whether the finite width 



of the internal propagators alone can provide the necessary cut-off fi, in particular when 
the singular contributions are evaluated on the corrected quasi-particle mass-shells, i.e. 
including damping. This has been cleared up recently in Ref. [J25|J . Our approach, however, 
was a strictly perturbative one. Through eq. (|35|) it requires to evaluate the corrections 
at the location of the (real) leading-order position, which, at this level, is in accordance 
with the findings of Ref. f25| . Our results and conclusions should therefore be quite inde- 
pendent of the details of the actual higher-order effects, as long as they are indeed able 
to cut off the quasi-particle mass-shell singularities. 

Acknowledgments: F. F. is supported by Deutsche Forschungsgemeinschaft (DFG). 



Appendix 



In the following we describe some of the steps which are encountered when deriving 
the main results of sects. 2 and 3 in a more pedestrian way. 
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The complete next-to-leading order correction to the longitudinal gluon polarization 
function 5Pl£, on the longitudinal mass-shell, has recently been written down in full detail 
in eq. (4.5) of Ref. [H5], albeit still at a purely algebraic level. To cover also the transverse 
case, we generalize the result obtained in and give an expression for Sn^ which is 
valid when used under the trace either with the matrix B^ v introduced in sect. 3 and 
taken at the longitudinal mass-shell or with Au V and taken at the transverse mass-shell: 



MP" = g 2 N J2 (eg" + Aj A t <% + Aj A t <% + A~A t c; 
with the coefficient matrices & v 



p2 K 2 



2P 2 g^ - 4pA t P t 



p2 

■p/xOO *yu00 _ r r- 

2p 2 k 2 2p 2 1 



p2 K 2 R 2 

k 2 p 



p 2 k 2 

p2 R 2 



P 2 S 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



2p 2 k 2 



K 2 
k 2 



p2 K 2 

Q-P,k 



P^ + RtR h 



2p 2 k 2 

q - p, A 



(A.5) 



The notation is that of @ : K = Q - P, k = q - p, A = 1/P 2 , Aj = 1/(P 2 - Pi(P)) 
(z = i, t), 5i = P 2 - ili(P) = A" 1 , U = (1,0) . An index minus refers to the shift 
P — > K in the corresponding quantity. The four vector R has no zeroth component : 
R» = (<5 t -^P 7P 2 )(P [P-P^). The arguments of the 4-leg vertex T are Q, -Q, -P, P , 
those of all 3-leg T's are Q, —K, —P . Finally, the unadorned sum symbol is short-hand 



for E P Jd 3 p/(2 



7T 



The above result ( |A.lj ) has a convergent sum over P , and its p-integration is restricted 
to soft values automatically due to the subtraction of the hard contribution. The calcu- 
lation was done in general covariant gauge with gauge parameter a, which dropped out 
algebraically. 

As in sect. 3, we concentrate on contributions involving a transverse propagator, 

£A t Ar c (P ,p) (A.6) 

Using the spectral representation 

Pt(x,p) 



dxx 



P 2 - x 2 



(A.7) 
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we may write 

d 3 p ™ 1 ~ 2 



j^J_Jx-p t (x lP )Y, 1 ^—^A t (Q -P ,c l -p)c(P ,p) . (A.8) 



Po 



Inspecting now the integration region of small p, one realizes that the weight of the 
transversal (but not the longitudinal) density (times 1/x) is concentrated at x = (see 
e.g. eq. (B.13) of Ref. pi): 



-Pt(x,p) -> \s(x) (p 2 <m 2 ) . (A.9) 

x p z 



Inserting ( |A.9| ) in (|A.8| ) and performing the ^-integration, one finds that the sum over P 



reduces to the term P = provided c(P , p) has no poles at P = 0, but this is easily 
excluded by inspection. Infrared singularities can now occur when A~ or c(0, p) diverges 
for p — > 0. By inspection one finds that c(0, p) is regular, so all singularities are due to 
the mass-shell singularities when approaching mass-shell i. The first term of the Taylor 
series for c(0, p) is responsible for the dominant contribution, which thus involves only 
the paradigmatic term Si studied in sect. 2 times functions of the external momentum. 



This way eqs. flT3) and (GDI) are readily reproduced 
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